Stabilization of the Cart–Inverted-Pendulum System Using State-Feedback Pole-Independent MPC Controllers

In this paper, a pole-independent, single-input, multi-output explicit linear MPC controller is proposed to stabilize the fourth-order cart–inverted-pendulum system around the desired equilibrium points. To circumvent an obvious stability problem, a generalized prediction model is proposed that yields an MPC controller with four tuning parameters. The first two parameters, namely the horizon time and the relative cart–pendulum weight factor, are automatically adjusted to ensure a priori prescribed system gain margin and fast pendulum response while the remaining two parameters, namely the pendulum and cart velocity weight factors, are maintained as free tuning parameters. The comparison of the proposed method with some optimal control methods in the absence of disturbance input shows an obvious advantage in the average peak efficiency in favor of the proposed SIMO MPC controller at the price of slightly reduced speed efficiency. Additionally, none of the compared controllers can achieve a system gain margin greater than 1.63, while the proposed one can go beyond that limit at the price of additional degradation in the speed efficiency.


Introduction
The cart-inverted pendulum (CIP) system that belongs to the class of fast singleinput, multiple-output (SIMO), under-actuated systems and satisfies a set of complicated characteristics, such as fourth-order highly nonlinear dynamics, open-loop instability, state coupling, and non-minimum-phase (NMP) behavior, provides many challenging problems to standard and modern control techniques [1]. In the context of the CIP system stabilization, moving the cart from an initial position to a final destination while keeping the pendulum erected in the upright position has been extensively studied in the past, and many output-feedback and (static and dynamic) state-feedback control techniques have been developed to solve it. However, solving this task efficiently in the framework of linear static state-feedback control (SFC) to ensure prescribed system gain margin in addition to good time response behavior with pole-independent parameter tuning is a subject that still needs further investigation.
There are different types of control methods that have been applied to the inverted pendulum systems [1], including model predictive control (MPC) and non-MPC methods. With regard to the complicated characteristics of the inverted pendulum plants, the needed controlled system performance, and the limited control input effort resource, time-domain optimization techniques, such as the MPC [2][3][4][5][6], seem to be one of the most convenient ways to tackle the above control problem, especially when state and control input constraints are considered. The key feature of the MPC method is based on the following three successive steps [3]: (i) the explicit use of a model and system measurements to predict proposed for automatic parameter tuning [24]. Notice that with a reduced set of parameters, the tuning difficulty becomes less problematic. Therefore, if a controller with few tuning parameters shows good performance for the considered problem, this can be seen as a huge practical advantage.
As can be seen, there are several interesting attempts to design linear static SFC for the CIP system stabilization in the form of MPC or non-MPC control schemes. MPC control schemes are more significant, because they can be considered optimal for a specified cost function. Notice that an ECS can be considered as an MPC method if there is a correspondence between the ECS gains and MPC parameters. The comparison, conducted in [25], between the MPC and LQR has shown that the MPC method is more suitable for the trajectory tracking task and smoothing in the control input, while the LQR is more convenient for fixed-value control and disturbance rejection, but it may generate adverse and rapid changes in the control signal. However, for both approaches, MPC and non-MPC, the presence of real NMP zeros in the cart part of the fourth-order linearized CIP transfer function limits the robustness performance and prevents the achievement of monotonic cart step responses. In order to obtain the best possible performance, the optimal choice of the SFC gains needs to be addressed. An important contribution in this line for pole-dependent non-MPC SFC controllers is attributed to the authors of [17], who proposed an analytical formula for optimal tuning of the SFC gains for the CIP system. In the derivation of the formula, the authors promote a priori a coincident-pole structure, which has a single tuning parameter (see Appendix C) for the closed-loop poles before maximizing the worst gain margin associated with the CIP output signals. In doing so, it is clear that the adopted configuration will impose an upper limit for the achievable system gain margin and prevent the controller from exploiting other possible pole configurations that may be more helpful in specifying a priori a prescribed gain margin and in reducing the impact of the closedloop CIP zeros on its performance. In addition to the CIP stabilization problem, designing controllers that achieve non-overshooting/undershooting for all-pole systems (i.e., systems only having poles in their transfer functions) or minimum overshooting/undershooting for non-all-pole systems (i.e., systems having poles and zeroes in their transfer functions) have received considerable attention [26][27][28][29][30]. In [26], for example, the authors studied the overshoot of an all-pole fourth-order system with respect to the variation of pole locations, where the poles are parameterized with two damping ratios and two undamped natural frequencies. One of the main obtained results states that the system step input overshoot remains unchanged if the ratio of their two natural frequencies is kept constant. The other main result states that the overshoot of the considered system does not have monotonicity with respect to each damping ratio or their sum. The above statements no longer hold for our CIP system, which is a non-all-pole system. In this case, and due to the presence of a single real NMP zero, achieving the cart step response monotonicity with linear SFC is almost impossible [29,30], and it appears natural, as a possible objective, to turn toward the design of controllers that achieve as little undershoot/overshoot as possible while maintaining good robustness stability.
On the other hand, addressing the SFC gains tuning problem to achieve monotonic step responses with pole-independent tuning methods can be performed for all-pole systems using the well-known coefficient diagram method (CDM) [27,28,31]. In this method, controllers are designed via the assignment of the so-called characteristic ratios and generalized time constant (GTC), which may have a strong physical relationship with the damping (i.e., overshoot) and speed of response of the closed-loop system, respectively. With non-allpole systems, the complexity of pole-zero interaction makes the standard CDM no longer valid. However, in the case of a non-all-pole system with one pair of jω-axis zeroes, the authors of [28] showed the possibility to obtain monotonic step response if the GTC is kept above a certain lower bound. For the CIP system, the cart system can be decomposed into the difference between a non-all-pole system with one pair of jω-axis zeroes and an all-pole system. Since it is impossible to ensure the cart monotonic step response, one can suggest as another possible objective to ensure the monotonicity of the non-all-pole system Sensors 2022, 22, 243 4 of 28 in the context of Manabe form, as is done in [28], without considering the impact of the remaining all-pole system (see Appendix A). In such a case, the robustness stability issue must be addressed, as well as the overshoot and undershoot response requirements (see Appendix A).
To the authors' knowledge, the attempts to meet a prescribed system gain margin, short CIP settling time with insignificant overshoot and undershoot, and reduced control effort with pole-independent parameter tuning have not been considered before in the context of linear static SIMO MPC SFC control schemes. The proposed MPC controller is designed with a generalized prediction model to circumvent an obvious stability problem [32] and a quadratic cost function with four control parameters. The first two parameters, namely the horizon time and the relative cart-pendulum weight factor, are automatically adjusted to ensure a priori prescribed system gain margin and fast pendulum response, while the remaining two parameters, namely the pendulum and cart velocity weight factors, are maintained as free tuning parameters to tackle the damping and the reduced control effort problem. In contrast to the existing state-feedback methods to stabilize the CIP system, this work has the following main distinguishing features:

•
From the theoretical point of view, a new robust pole-independent SIMO MPC controller with only two adjusted control parameters is proposed to solve the fourth-order CIP stabilization problem under full state availability, known CIP parameters, and a pendulum mass that is negligible in comparison to the cart mass. It is shown in this paper how to constraint the SFC controllers to ensure a priori a prescribed CIP system gain margin and how to constrain the MPC controller to have a fast pendulum response in comparison to the cart response, i.e., satisfying a two-time-scale structure in which the closed-loop pendulum subsystem responds faster than the closed-loop cart subsystem [33]. This contribution leads effectively to the reduction in the number of free tuned parameters from four to two. Another contribution, mainly inspired by the work of [17], is to derive optimal controllers that maximize the system gain margin for the standard CDM and cascade methods (see Appendices A-C). Finally, in order to obtain the optimal MPC weight factors and to allow performance comparison, the indices of speed and average peak efficiencies are introduced for characterizing the closed-loop CIP transient responses.

•
From a practical point of view, for a prescribed system gain margin, the impact of the velocity cost function weight factors on the closed-loop system transient performance in the context of the proposed MPC method can be evaluated off-line, and the obtained trends can be easily clarified using two-dimensional graphical and contour plot representations. Some useful guidelines for rapid weighting factor adjustment are developed for the proposed method in the presence (or absence) of disturbance input. Additionally, standard and advanced numerical tuning algorithms can benefit from the obtained reduced two-dimensional space search to achieve global optimality for a given criterion. Such a situation helps in rapidly obtaining a solution to the parameter tuning problem and also helps in checking easily its optimality.
The paper is organized as follows. Section 2 defines the CIP models and states the problem under consideration. Section 3 deals with the SIMO MPC controller design. Section 4 provides simulations results, and finally, Section 5 concludes this paper.

Nonlinear Inverted Pendulum Dynamics
The CIP system consists of a cart and a rigid rod pendulum with a pivot mounted on the top of the cart, as shown in Figure 1. Under the action of the horizontal force that is regarded as the control input u(t), the cart moves left or right on a one-dimensional bounded track, whereas the pendulum swings in the vertical plane determined by the track. It is assumed that no friction exists in the system between the cart and the track or between the cart and the pendulum. The cart is characterized by a mass M, and the pendulum is track. It is assumed that no friction exists in the system between the cart and the track or between the cart and the pendulum. The cart is characterized by a mass M , and the pendulum is characterized by a length L of the mass-less rod and a point mass m located at the free end of the pendulum rod. To describe the dynamics, a state vector is used, where the state signals 1 x , 2 x , 3 x , and 4 x are the pendulum angle (measured clockwise w.r.t the upward vertical), the angular velocity of the pendulum, the cart position, and the cart velocity, respectively. Using Newton's laws, it is possible to derive and describe the dynamics of the CIP system explicitly by the following fourth-order SIMO under-actuated and highly nonlinear dynamical system (see Appendix D): ( ) (1) Notice that, in the above model, the pendulum and cart accelerations are driven with opposite direction input control forces over the upper half-plane of the pendulum angle. This makes our stabilization task more challenging.

Generalized Linear Prediction Model
Under the three assumptions taken from [12], the assumption of the light mass at the end of the rod, and the control input domination assumption the nonlinear dynamical model (1) can be approximated by the following linear model: with Figure 1. The inverted pendulum system.
To describe the dynamics, a state vector x = (x 1 , x 2 , x 3 , x 4 ) T is used, where the state signals x 1 , x 2 , x 3 , and x 4 are the pendulum angle (measured clockwise w.r.t the upward vertical), the angular velocity of the pendulum, the cart position, and the cart velocity, respectively. Using Newton's laws, it is possible to derive and describe the dynamics of the CIP system explicitly by the following fourth-order SIMO under-actuated and highly nonlinear dynamical system (see Appendix D): (1) Notice that, in the above model, the pendulum and cart accelerations are driven with opposite direction input control forces over the upper half-plane of the pendulum angle. This makes our stabilization task more challenging.

Generalized Linear Prediction Model
Under the three assumptions taken from [12], the assumption of the light mass at the end of the rod, and the control input domination assumption the nonlinear dynamical model (1) can be approximated by the following linear model: with Notice that the parameters a 1 , b 1 , and b 2 are positive. It will be shown later in the paper that all the assumptions given in (2) are indeed satisfied. From (3), one can deduce the following unstable transfer functions: Let us now define the following generalized linear CIP model with k ∈ {−1, +1}, which can be seen as a generalization or modification of model (3): where the pendulum subsystem model remains unchanged when compared to the original model (3) while the cart subsystem model is substituted by a fictitious one. This last subsystem, with the new states y 3 and y 4 , is driven at each instant t with the control input ku. Using the above-adopted model (6), and assuming y 3 = x 3 and y 4 = x 4 at the instant t, we define the following pendulum and cart generalized prediction models: where h is a positive constant horizon time. If u(t) = u is constant on the interval [t, t + h], then by solving (6) and using the last equality of (4), we obtain (7) with the following parameters: Notice that the parameters A 1 , A 2 , B 1 , B 2 , E 1 , E 2 , E 3 , and E 4 are positive.

Problem Statement
For the CIP system stabilization, we are interested in using the following static statefeedback control law (SFC): where N = (N 1 , N 2 , N 3 , N 4 ) is the control gain vector to be determined, and x 3d is the constant cart reference position. The main result of this paper is how to simplify the control tuning based on the good robustness properties of the control system and the optimized two-time controller structure. The controllability of the linear system (3) ensures the existence of a set of SFC having the form (9) that can achieve the stabilization of the CIP system in the vicinity of the unstable equilibrium point. Combining (5) and the Laplace transform of (9) yields: The obtained linearized closed-loop CIP system is of the fourth order and has the basic configuration shown in Figure 2. Obviously, each controller belonging to the class (9) can be interpreted as the combination of two PD controllers (one for the pendulum subsystem and the other for the cart subsystem).
The obtained linearized closed-loop CIP system is of the fourth order and has the basic configuration shown in Figure 2. Obviously, each controller belonging to the class (9) can be interpreted as the combination of two PD controllers (one for the pendulum subsystem and the other for the cart subsystem).

Remark 1.
The fourth-order transfer functions (10) must at least be stable. The Routh necessary and sufficient stability conditions for such a system are given in (11). From these conditions and the positivity of the parameters (4), it follows that the controller gains ( )   , provokes the appearance of an undesirable initial undershoot in the cart step response. The amplitude of this undershoot grows to infinity if the settling time is reduced to 0 [30]. This behavior therefore places a limitation on the cart speed of response. Overshoot in the

Remark 1.
The fourth-order transfer functions (10) must at least be stable. The Routh necessary and sufficient stability conditions for such a system are given in (11). From these conditions and the positivity of the parameters (4), it follows that the controller gains (N 1 , N 2 , N 3 , N 4 ) must be positive.
Remark 2. The cart transfer function F 3 (s) has two real opposite zeroes, z 1,2 = ±a 1/2 1 . The single real NMP zero, z 2 = +a 1/2 1 , provokes the appearance of an undesirable initial undershoot in the cart step response. The amplitude of this undershoot grows to infinity if the settling time is reduced to 0 [30]. This behavior therefore places a limitation on the cart speed of response. Overshoot in the above response is another undesirable effect that may be reduced with the undershoot phenomenon if an appropriate selection of SFC controller gains is conducted. Remark 3. From (10), it can be seen that the pendulum and cart transfer functions F 1 (s) and F 3 (s) have fixed zeroes and adjustable poles. Thus, whatever the method used to determine the SFC gains, it is always considered as a pole placement method. In the sequel, we shall refer to the SFC control method (9) as pole-dependent and pole-independent if its gain tuning method specifies a priori and a posteriori the closed-loop pole locations of (10), respectively.
In this paper, we investigate a subclass of SFC controllers (9), where the designed controller carries all the key features of SIMO MPC controllers and can be implemented in an ECS fashion. As it is depicted in Figure 3, the considered controller receives as its inputs the cart reference position, x 3d , and the state vector, x = (x 1 , x 2 , x 3 , x 4 ), and produces as its single output the control signal, v, which is fed directly to the control input u of the CIP system. In our investigation, we need to consider the following assumption. have fixed zeroes and adjustable poles. Thus, whatever the method used to determine the SFC gains, it is always considered as a pole placement method. In the sequel, we shall refer to the SFC control method (9) as pole-dependent and pole-independent if its gain tuning method specifies a priori and a posteriori the closed-loop pole locations of (10), respectively.
In this paper, we investigate a subclass of SFC controllers (9), where the designed controller carries all the key features of SIMO MPC controllers and can be implemented in an ECS fashion. As it is depicted in Figure 3, the considered controller receives as its inputs the cart reference position, d x 3 , and the state vector, ( ) , and produces as its single output the control signal, v , which is fed directly to the control input u of the CIP system. In our investigation, we need to consider the following assumption.
Assumption 1: • the state vector ( ) are known constants; and • the set of hypotheses in Equation (2) is satisfied. Under Assumption 1, the CIP stabilization problem may be formulated in two steps as follows. In the first step, restrict the class of SFC controllers so as to satisfy the following property: • Property 1: The SFC controller (9) ensures robust stability with a prescribed closed- Assumption 1: the parameters (M, m, L, x 3d ) are known constants; and • the set of hypotheses in Equation (2) is satisfied.
Under Assumption 1, the CIP stabilization problem may be formulated in two steps as follows. In the first step, restrict the class of SFC controllers so as to satisfy the following property: • Property 1: The SFC controller (9) ensures robust stability with a prescribed closed-loop gain margin, GM min .
In the second step, and under the above gain margin constraint and the generalized prediction model (7) and (8), design a pole-independent SIMO MPC controller such that the following properties are satisfied: • Property 2: The closed-loop CIP system satisfies a two-time-scale structure in which the closed-loop pendulum subsystem responds faster than the closed-loop cart subsystem [33]. • Property 3: In the absence of disturbance input, the controller ensures low cart settling time T > 0 without excessive peaking (undershoot and overshoot) phenomenon in the CIP input and output responses.

SIMO MPC Controller Design
To design robust SIMO MPC controllers, we first developed Equation (23), associated with the SFC gain margin property, using the stability conditions (11) and the prescribed gain margin, GM min . Then, Property 1 holds when the above-mentioned equation is satisfied. Next, based on the generalized prediction model (7)(8) and the quadratic cost function (24) that is characterized by a set of MPC parameters ϕ = (h, ρ 1 , r, ρ 2 ), i.e., the horizon time and the positive weight factors to be defined later, we establish the explicit relationship (26) between ϕ and N. Finally, we propose a pole-independent tuning method that relies primarily on choosing the parameter set ϕ of the cost function as a starting point of the optimal control. This method leads to determining the pole locations a posteriori, i.e., at its final stage. Notice that the controller design developed hereafter uses the linearized CIP model (3) in the neighborhood of the equilibrium point, where it is assumed to be justified at least from the approximation point of view. The validity of the adopted design methods is demonstrated from the control point of view in Section 4, where the nonlinear model (1) is used instead of its linearized version.

Closed-Loop SGM Constraint
Let us consider the following change of variables: where α 12 > 0 and α 34 > 0. Introducing (12) in the last condition of (11) gives: This means that there exists a critical gain N 4 that depends on the fixed value of N 2 . The critical positive gain b 2 N 4i > 0 at which the system becomes marginally stable is: Sensors 2022, 22, 243 9 of 28 From (14), we define the cart-loop gain margin as follows: On the other hand, the critical positive gains b 1 N 2i α > 0 and y i = b 1 N 2i α − a 1 at which the system becomes marginally stable is deduced from (13) as the positive solutions of the following system of equations: with Using the last equality gb 1 = a 1 b 2 of (4), (17) reduces to: From (18), we define the pendulum-loop gain margin as follows: Taking into account the worst case, the system gain margin index can be defined as follows [17]: Since the gain margins, GM C and GM P , must be greater than 1, from (16) and (19), we deduce, respectively: Regarding (21), the system gain margin (20) reduces to: The following remark states some of the obvious properties of the closed-loop system gain margin.

Remark 4.
Consider the fourth-order closed-loop system in (10) with the change of variables (12) and the system gain margin (22). Then, the following properties hold: (i) with fixed parameters α and b 1 N 2 , GM decreases with the increase in b 2 N 4 ; (ii) with fixed parameters α and b 2 N 4 , GM increases with the increase in b 1 N 2 ; and (iii) with fixed parameters b 1 N 2 and b 2 N 4 , GM remains unchanged for any value γ if N 1 and N 3 are increased by the amounts N 2 γ and N 4 γ, respectively.

Remark 5.
Given the value of GM min and using Equation (22), Property 1 can be satisfied under the following constraint: where α is defined in (12).

Deriving the SIMO MPC Control Law
To derive the SIMO MPC control law, we first adopt the generalized linear prediction model (7) and (8) as a prediction model for the CIP system. Then, we consider the following quadratic cost function: where ϕ = (h, ρ 1 , r, ρ 2 ) ∈ +4 is the set of positive MPC parameters, e 1 (t + h) and e 2 (t + h) are the predicted pendulum angle and velocity errors, while e 3 (t + h) and e 4 (t + h) are the predicted cart position and velocity errors. Now, given the cart destination x 3d and the prediction models (7), we obtain the optimal input v 0 , which minimizes the value of the quadratic cost function (24). Recall that the references associated with x 1 , x 2 , and x 4 are assumed to be zero. Substituting (7) into (24) and setting the gradient of J MPC (t, h) with respect to v to zero yields: with the SIMO MPC controller gains given below: Notice that the positivity of (ρ 1 , r, rρ 2 ) together with the positivity of (8) imply the positivity of N 1 and N 2 in addition to the fact that the gains N 3 and N 4 have an undefined sign that follows the one associated with the value of k. Remark 6. Since the signs of N 3 and N 4 are identical to the sign of k, the designed SIMO MPC controller that is defined by (23) and (24) does not ensure the stability condition (11) when using the trivial prediction model, i.e., the generalized prediction model (10) with k = +1. To circumvent such a drawback, we put k = −1 in the above model. This modification is equivalent to the substitution of the original cart subsystem model by another cart subsystem model that is driven at each instant and from the same state point with an opposite control input sign. The modification of the cost function while keeping the trivial prediction model is another option to solve the encountered stability problem. However, this option is more complicated and is not considered in this paper.

Time-Scale Structure Constraint
Here, we are interested in designing a SIMO MPC controller that ensures Property 2, i.e., a two-time-scale structure, without using a priori targeted closed-loop pole locations. To this end, the CIP control system of Figure 2 is transformed to the configuration of Figure 4, where a virtual reference x 1d and a new gain K are introduced for the closed-loop pendulum subsystem. The pendulum transfer function to the considered reference and the value of K that leads to a unity static gain are given by: where PL τ is a lower bound for the pendulum GTC.  The GTC associated with (10) and (27) are given, respectively, by: where τ PL is a lower bound for the pendulum GTC. For the sake of developing a general tuning method that is valid for a large class of CIP systems, let us define the following normalized parameters: where ϕ 0 is the set of normalized positive MPC cost function parameters. Substituting (29) into (28) and (26), we obtain: with Now, combining (30) and (32) yields: In order to speed the pendulum response, we plan to maximize (33) by a proper choice of the horizon time. More precisely, for a given positive value ρ 10 , setting the gradient of r 12 with respect to h 0 to zero yields: From (34) and Figure 5, it is clear that the parameter ρ 10 increases with the increase in h 0 and satisfies 0 < ρ 10 < 1. Concerning the amount r 12 , it is obvious that this parameter decreases toward 1 with the increase in h 0 . ρ remains to be determined. To this end, it is reasonable to satisfy Property 2 by assuming that the GTC of the closed-loop cart subsystem is lower than the GTC of the pendulum subsystem. According to (30) and (32) and the fact that 0 0 > r , this leads us to impose the constraint: , the cart settling time at 5%. To derive general comments on the behavior of the CIP system, it is useful to introduce  (23) gives Substituting the expressions of N 40 and D 0 from (32) in (35) and solving for r 0 gives: Since h 0 and ρ 10 are linked by the relationship (34) and the fact that r 0 and (h 0 , ρ 10 , ρ 20 ) are linked by the relationship (36), we only need to specify two parameters among the elements of the set ϕ 0 to achieve the controller design. In our case, we shall fine-tune the weight factor vector ρ 0 = (ρ 10 , ρ 20 ) inside an appropriate domain. We have already shown that 0 < ρ 10 < 1, and the interval associated with ρ 20 remains to be determined. To this end, it is reasonable to satisfy Property 2 by assuming that the GTC of the closed-loop cart subsystem is lower than the GTC of the pendulum subsystem. According to (30) and (32) and the fact that r 0 > 0, this leads us to impose the constraint: from which we obtain: A lower bound ρ 20min , which is shown in Figure 5 as a function of h 0 , is then imposed to ρ 20 to ensure Property 2. Now, let us define a set of five indices to evaluate the CIP transient response performance: P x1 (ρ 0 ), the maximum absolute value of the pendulum angle response; V x3 (ρ 0 ), the cart response overshoot; D x3 (ρ 0 ), the cart response undershoot; P u (ρ 0 ), the maximum absolute value of the control input signal; and t cs (ρ 0 ), the cart settling time at 5%. To derive general comments on the behavior of the CIP system, it is useful to introduce the change variable s = a 1/2 1 z and the gains (31) in (10) to obtain the following transfer functions: With the above processing, all the considered CIP transient response characteristics remain unchanged, except the one associated with the modified (normalized) cart settling time t csn , which is now linked to the original cart settling time t cs by the relationship: Regarding the transfer functions (39), it appears naturally useful to normalize the indices associated with the pendulum angle response and control input signal as follows: Now, we are ready to summarize in Algorithm 1 how to design the proposed SIMO controller and how to obtain the closed-loop linearized CIP performance when the system gain margin GM min and the weight vector ρ 0 = (ρ 10 , ρ 20 ) are chosen. Algorithm 1: Controller design and closed-loop linearized CIP system evaluation.
Evaluate the indices (P x1n , V x3n , D x3n , P un , t csn ) from the obtained step responses; 7. End.

Parameter Tuning
In this section, parameter tuning for the weight vector ρ 0 is performed to ensure the validity of Property 3. To this end, the peaking and speed constraints need to be considered simultaneously. To solve this problem, a time-domain optimization technique that involves tackling two issues, namely the choice of the performance criterion and the mean to optimize it, is proposed hereafter. For the first issue, we have found it useful, regarding the desirable Property 3, to define two proposed sound and scaled-free indices. These indices are the speed efficiency S E (ρ 0 ) and the average peak efficiency S E (ρ 0 ) and are defined for a given standard (reference) state-feedback method (SSF) and our control method as follows: where P 1 (ρ 0 ) = P x1 (ρ 0 ), P 2 (ρ 0 ) = P x3 (ρ 0 ), P 3 (ρ 0 ) = P u (ρ 0 ), and P 4 (ρ 0 ) = t cs (ρ 0 ). The index P i,SSF has the same interpretation as the index P i (ρ 0 ), with the proposed control method replaced by the SSF one. Obviously, an efficiency index greater than 50% indicates that the proposed control method outperforms the SSF method; otherwise, degradation in the performance of our method in comparison to the SSF one is noted. Then, taking into account the worst efficiency case, we may formulate the setting problem as a maximin optimization model as follows: where λ = 0 if speed efficiency is the only concern, and λ = 1 if peaking and speed efficiency are both considered in the optimization. It should be noted that (43) is a continuous nonlinear optimization model with a highly nonlinear (possibly discontinuous) objective function J 1 (ρ 0 ), and it is difficult to know a priori whether such a function is unimodal or multimodal before starting the optimization. To avoid erroneous solutions, problem (43) has to be solved to global optimality in the considered parameter space domain. Then, to solve the second issue, we opt for global optimization using an exhaustive search over the domain that is defined by 0 ≤ ρ 10 ≤ ρ 10,max and 0 ≤ ρ 20 ≤ ρ 20,max . This is interesting, since there are only two tuning parameters.

Numerical Simulations
Numerical simulations are conducted in several separate sections to show the potential of the proposed pole-independent SIMO MPC controller (SIM) and its advantage in comparison to the pole-independent standard CDM controller (CDM), the pole-dependent CASC MPC controller (CAS), and the coincident pole controller (CPP). The definitions of the CDM, CAS, and CPP controllers are given in Appendices A-C, respectively. These controllers are chosen for their efficiencies and for the reduced number of adjusted parameters, which leads to yield, without significant difficulty, a guarantee of the overall best performance for each tuned controller. Based on the linearized CIP system, we begin the numerical simulation in Section 4.1 with the establishment of some guidelines that help in choosing the SIM weight factors for a prescribed system gain margin in the absence of disturbance input. These guidelines are quite general, since they can be applied to any linear CIP model of the form (3), as long as Assumption 1 stays satisfied. Additionally, they are very useful since they inform us how to tune the SIM weight factors to get some desirable transient closed-loop linearized CIP system performance. Next, to make a fair comparison, we tune, for a given physical CIP system, the best possible CDM, CAS, and CPP controller so as to obtain a maximum system gain margin for each tuned controller. For the SIM controller, the system gain margin is set a priori as the best one so far obtained by the above controllers, and the SIM weight factors are tuned according to (43). In Sections 4.3 and 4.4, we compare the obtained controllers on the nonlinear CIP system in the absence and presence of disturbance input.

Guidelines for Weighting Factor Adjustment
To obtain some insights into how to choose the weight vector ρ 0 = (ρ 10 , ρ 20 ) for different system gain margin values in the absence of disturbance input, we conducted 50,000 simulations with Algorithm 1 (obtained by using a coarse grid discretization with 50 points for ρ 10 in the interval 0 < ρ 10 < 0.5, 50 points for ρ 20 in the interval 0 < ρ 20 < 5, and 20 points for GM min in the interval 1 < GM min < 3). Figures 6 and 7 show the obtained contour plots for the considered transient response characteristics for GM min = 1.63. To generate the contours, we constructed a grid interpolant using the Matlab function "griddedInterpolant" with the "pchip" option and interpolate the considered index with 10 −3 spacing. To avoid dummy solutions, only the contours that have more than 10 3 points were retained. the cart speed of response, there are two separate regions, i.e., the upper region UR and the lower region LR , where the settling time can be reduced below the value 20. These regions are located above the red line, where the increase in 20 ρ while maintaining 10 ρ constant reduces the peaking phenomenon. Table 1 shows the associated region frontier locations and transient performance intervals. Obviously, choosing 0 ρ in the UR appears to be more recommended than choosing it in the LR .    When the cart speed is the only concern, i.e., when using the proposed tuning method (43) with 0 = λ , there is no need to select an SSF method, and the best-achieved cart settling time, (44) Figure 8 shows the above indices together with the cart overshoot and undershoot as a function of the prescribed system gain margin. It is clearly seen that the best-achieved cart settling time is obtained with 2 . 1 min ≈ GM ; additionally, the increase in min GM above this value leads to an increase in the optimal cart settling time and a reduction in the cart undershoot. The cart overshoot remains between 4 and 5%, while the optimal weight factor b , 10 ρ remains, in all studied cases, under 0.12. As we shall see in the next section, the considered controllers for the comparison task do not exceed a system gain margin of 63 . 1 , while our SIM controller can go beyond this limit, as is indicated in Figure 8. From Figures 6 and 7, one can observe that the choice of the weight factors cannot be performed independently to address the peaking phenomenon and the cart speed of response issues simultaneously. Additionally, a strong correlation between the contours of P x1n and D x3 is noticed. For the peaking phenomenon, it is clearly seen that the increase in ρ 20 while maintaining ρ 10 constant reduces without ambiguity the cart overshoot V x3 (ρ 0 ). Ensuring a non-overshooting behavior for the cart subsystem thus appears possible with a high enough value of ρ 20 . For a constant ρ 10 , the peaking phenomenon indices, namely P x1n , D x3 , and P un , are only reduced when ρ 0 is located above their associated red line frontiers, as shown in Figures 6 and 7. The evaluation of these frontiers is performed without considering the interpolation processing. Ensuring a reduced undershooting behavior for the cart subsystem appears to be possible with the increase in ρ 20 . In addition to that, one can also clearly observe that the increase in ρ 10 while maintaining ρ 20 constant contribute, without ambiguity, to reduce the control input effort P un and to increase the cart overshoot V x3 (ρ 0 ), as long as the value of ρ 10 still far enough from the frontier that it defines the validity of the MPC controller (see step 3 of Algorithm 1). For a constant ρ 20 , the increase in ρ 10 may lead to an increase in P x1n and D x3 . Concerning the cart speed of response, there are two separate regions, i.e., the upper region UR and the lower region LR, where the settling time can be reduced below the value 20. These regions are located above the red line, where the increase in ρ 20 while maintaining ρ 10 constant reduces the peaking phenomenon. Table 1 shows the associated region frontier locations and transient performance intervals. Obviously, choosing ρ 0 in the UR appears to be more recommended than choosing it in the LR. When the cart speed is the only concern, i.e., when using the proposed tuning method (43) with λ = 0, there is no need to select an SSF method, and the best-achieved cart settling time, t csn,b , and its associated weight factor, ρ 10,b , are defined as follows: t csn (ρ 10 , ρ 20 , GM min ) .
(44) Figure 8 shows the above indices together with the cart overshoot and undershoot as a function of the prescribed system gain margin. It is clearly seen that the best-achieved cart settling time is obtained with GM min ≈ 1.2; additionally, the increase in GM min above this value leads to an increase in the optimal cart settling time and a reduction in the cart undershoot. The cart overshoot remains between 4 and 5%, while the optimal weight factor ρ 10,b remains, in all studied cases, under 0.12. As we shall see in the next section, the considered controllers for the comparison task do not exceed a system gain margin of 1.63, while our SIM controller can go beyond this limit, as is indicated in Figure 8.

Disturbance-Free Parameter Tuning
Now, let us consider the nonlinear CIP system (1) with a set of physical parameters  [23]. Since the behavior of the linearized CIP system is considered quite similar to the behavior of the nonlinear CIP system in the vicinity of the equilibrium point, the tuning of the controllers is thus based only on the linearized system. Figures 9 and 10 show the evolution of the system gain margin together with the associated linearized cart transient performance for the considered controllers. In contrast to non-undershooting cart response behavior, ensuring a non-overshooting cart response behavior is a feasible objective for all controllers. Reducing the cart response undershoot can be undertaken at the price of an increase in the cart settling time and/or deterioration in the system gain margin

Disturbance-Free Parameter Tuning
Now, let us consider the nonlinear CIP system (1) with a set of physical parameters M = 2.4 kg, m = 0.23 kg, L = 0.36 m, g = 9.81 m/s 2 , and a cart track length limited between ±0.5 m [23]. Since the behavior of the linearized CIP system is considered quite similar to the behavior of the nonlinear CIP system in the vicinity of the equilibrium point, the tuning of the controllers is thus based only on the linearized system. Figures 9 and 10 show the evolution of the system gain margin together with the associated linearized cart transient performance for the considered controllers. In contrast to non-undershooting cart response behavior, ensuring a non-overshooting cart response behavior is a feasible objective for all controllers. Reducing the cart response undershoot can be undertaken at the price of an increase in the cart settling time and/or deterioration in the system gain margin performance. Regarding the system gain margin trend, it is seen that the system gain margin is limited by about 1.4 for the CDM method and by about 1.6 for the CAS and CPP methods. For these methods, dependencies are typically observed between the system gain margin and the tuning parameters, which render them less flexible. Now, following [17], we retain the best tuning parameter for the CDM, CAS, and CPP as those that maximize the system gain margin. The gains of these controllers and the resulting linear CIP system transient performance are summarized in Tables 2 and 3, respectively. The transient performances are obtained for zero initial conditions and a cart step x 3d = 0.1.

Disturbance-Free Parameter Tuning
Now, let us consider the nonlinear CIP system (1) with a set of physical parameters , and a cart track length limited between m 5 . 0 ± [23]. Since the behavior of the linearized CIP system is considered quite similar to the behavior of the nonlinear CIP system in the vicinity of the equilibrium point, the tuning of the controllers is thus based only on the linearized system. Figures 9 and 10 show the evolution of the system gain margin together with the associated linearized cart transient performance for the considered controllers. In contrast to non-undershooting cart response behavior, ensuring a non-overshooting cart response behavior is a feasible objective for all controllers. Reducing the cart response undershoot can be undertaken at the price of an increase in the cart settling time and/or deterioration in the system gain margin performance. Regarding the system gain margin trend, it is seen that the system gain margin is limited by about 1.4 for the CDM method and by about 1.6 for the CAS and CPP methods. For these methods, dependencies are typically observed between the system gain margin and the tuning parameters, which render them less flexible. Now, following [17], we retain the best tuning parameter for the CDM, CAS, and CPP as those that maximize the system gain margin. The gains of these controllers and the resulting linear CIP system transient performance are summarized in Tables 2 and 3, respectively. The transient performances are obtained for zero initial conditions and a cart step (a) (b)   Table 3. Controller performance comparison using linearized CIP model.  For the proposed SIM method, one can impose a priori the system gain margin as a constraint to the controller. Using a gain margin of 1.63 and applying the proposed tuning method (43) with λ = 1, we obtain ρ 10,g = 0.09 and ρ 20,g = 1.5. The obtained controller gains and their associated performances are also given in Tables 2 and 3, respectively. Notice that the last assumption of Equation (2) is satisfied, since we have, for all controllers, N 1 >> mg. In comparison to the other control methods, the proposed SIM method has the best pendulum angle deviation, cart undershoot, and control input effort at the price of degradation in the cart overshoot and settling time. Using the CPP for the linearized CIP system as an SSF, we obtained the best G E at a price of a slight degradation in S E , as is indicated in Table 3. Figure 11 indicates that the obtained optimal SIM controller depends only on S E , since we always have S E < G E . In addition to that, Table 4 and Figure 11 show that the optimal SIM controller is characterized by a set of two complex conjugates poles: one of them has a real part that is equal to the CPP controller pole, and the other one has the highest real part, which somewhat explains why the SIM controller has the largest cart settling time. Finally, Figure 10 tells us that the proposed SIM method allows to further reduce the peaking phenomenon while maintaining the system gain margin constant by increasing ρ 20 above 1.5 and maintaining ρ 10 = 0.09, but this enhancement is followed by an increase in t cs . In the vicinity of ρ 20 = 2.5, the overshoot vanishes, and the cart settling time takes a value of about 4.8 s.

Performance Analysis without a Disturbance Input
To compare the performance of the above controllers, the reference cart position d x 3 is set to m 0 , and all the initial state values are set to zeros, except the initial cart position, which is set to ( ) The simulation results applied to the nonlinear CIP system (1), without considering disturbance input, are shown graphically in Figures 12 and  13 and numerically in Table 5. One can notice at this point that all the assumptions given

Performance Analysis without a Disturbance Input
To compare the performance of the above controllers, the reference cart position x 3d is set to 0 m, and all the initial state values are set to zeros, except the initial cart position, which is set to x 3 (0) = −0.1 m. The simulation results applied to the nonlinear CIP system (1), without considering disturbance input, are shown graphically in Figures 12 and 13 and numerically in Table 5. One can notice at this point that all the assumptions given by Equation (2) are satisfied. The strong similarity that exists between the performance of the closed-loop linearized CIP system in Table 3 and the performance of the closed-loop nonlinear CIP system in Table 5 for each controller confirms the potential of linear control theory in solving the considered stabilization problem. On the other hand, Figure 12 shows that the SIM method exhibits good performance in the control input effort, the pendulum angle, and the cart position undershoot at the price of relatively large cart position overshoot and settling time. Concerning Figure 13, the fact that the control input effort demand occurs essentially at the beginning of the CIP system stabilization for all the considered controllers is obviously distinguished. Equations (26) and (29). When the pendulum angle is too large, the linearized model loses its validity, and it is necessary to apply some nonlinear technique to bring the pendulum to the vicinity of the upright equilibrium state. Then, it is possible to switch to the proposed control.

Performance Analysis with a Disturbance Input
To check the sensitivity of the studied controllers to disturbance input, a Matlab-Sim- When the initial states begin from an equilibrium point and reach another one while the assumptions of Equation (2) are satisfied during the motion, the controlled nonlinear CIP system behaves similarly to the controlled linear CIP system. In such a situation, the performance of the controlled nonlinear CIP system can be naturally evaluated using the contour plots of Figures 6 and 7 for a gain margin of 1.63 and given MPC parameters. However, when the states are far away from the equilibrium point, the classical prediction model, defined by Equations (7) and (8) together with k = +1, gives only partially consistent estimation or no longer holds at all. In the former case, which is not covered by the numerical study, it may perhaps be interesting to design the SIMO MPC controller in such a way to increase enough ρ 0 = (ρ 10 , ρ 20 ) in the hope to simultaneously increase the pendulum gains N 1 and N 2 while reducing the cart gains N 3 and N 4 , as is suggested by Equations (26) and (29). When the pendulum angle is too large, the linearized model loses its validity, and it is necessary to apply some nonlinear technique to bring the pendulum to the vicinity of the upright equilibrium state. Then, it is possible to switch to the proposed control.

Performance Analysis with a Disturbance Input
To check the sensitivity of the studied controllers to disturbance input, a Matlab-Simulink band-limited white noise n(t) with sampling time 0.01 s, seed = 23,341, and power P = 10 −3 was added to the control input signal of the nonlinear CIP system over a time interval of 120 s [23]. This disturbance input, which has an amplitude between −1.17 and 1.25, is depicted in Figure 14. Table 6 summarizes the resulting maximum absolute values of the CIP input and output signals that are obtained from the analysis of the last 100 s, i.e., in the assumed steady-state region. Regarding the disturbance input amplitude that appears in Figure 14, the controller performance of Table 5, and the maximum absolute values of Table 6, all the considered controllers appear to exhibit useful noise rejection capabilities. The CAS method, with its highest gains (see Table 2) and high transient control input effort (see Table 5), exhibits the best noise rejection capabilities in the steady-state region, whereas the proposed SIM method, with its lowest gains (see Table 2) and lowest transient control input effort (see Table 5), exhibits the worst one in the same steady-state region. Figures 15 and 16 show graphically the obtained CIP output responses. Although the disturbance input level appears to dominate completely the amplitude of the SIM transient noise-free control effort, i.e., 0.54 N (see Table 5), the proposed SIM method still performs well even with a tuning method that does not take into account the presence of a disturbance input.
ble 2) and lowest transient control input effort (see Table 5), exhibits the worst one in the same steady-state region. Figures 15 and 16 show graphically the obtained CIP output responses. Although the disturbance input level appears to dominate completely the amplitude of the SIM transient noise-free control effort, i.e., N 54 . 0 (see Table 5), the proposed SIM method still performs well even with a tuning method that does not take into account the presence of a disturbance input.     ble 2) and lowest transient control input effort (see Table 5), exhibits the worst one in the same steady-state region. Figures 15 and 16 show graphically the obtained CIP output responses. Although the disturbance input level appears to dominate completely the amplitude of the SIM transient noise-free control effort, i.e., N 54 . 0 (see Table 5), the proposed SIM method still performs well even with a tuning method that does not take into account the presence of a disturbance input.    ). Figure 17 shows the obtained contour plots for the maximum absolute values of the input and outputs nonlinear CIP system. It Let us now study the impact of the weight vector ρ 0 = (ρ 10 , ρ 20 ) on the closed-loop nonlinear CIP steady-state performance when the adopted disturbance input is considered. To this end, we keep GM min = 1.63 and conduct 2500 simulations (obtained by using a coarse grid discretization with 50 points for ρ 10 in the interval 0 < ρ 10 < 0.5 and 50 points for ρ 20 in the interval 0 < ρ 20 < 5). Figure 17 shows the obtained contour plots for the maximum absolute values of the input and outputs nonlinear CIP system. It can be seen that the weight factor ρ 10 has a strong impact on σ x1 , σ x3 , and σ u , while ρ 20 has only a moderate impact on σ x3 . In addition, one can clearly see that the increase in ρ 10 reduces the CIP input index σ u , but as a side effect, it increases the CIP output indices σ x1 and σ x3 . Let us now study the impact of the weight vector ( ) 20 10 0 , ρ ρ ρ = on the closed-loop nonlinear CIP steady-state performance when the adopted disturbance input is considered. To this end, we keep 63 . 1 min = GM and conduct 2500 simulations (obtained by using a coarse grid discretization with 50 points for 10 ρ in the interval ). Figure 17 shows the obtained contour plots for the maximum absolute values of the input and outputs nonlinear CIP system. It can be seen that the weight factor 10 ρ has a strong impact on 1 Remark 7. This paper does not address the measurement noise issue explicitly. The experience with the cascade MPC controller in [17] suggests that low and moderate measurement noises are not problematic, and simple low-pass filtering can be used. This is also favored in [33]. However, high measurement noises present a challenge to the implementation of linear static state-feedback controllers. The problem lies in the need for differentiation of the pendulum angle and the cart position w.r.t. time. In such a situation, it is useful to consider sophisticated differentiations or other controller structures such as the output-feedback controller with or without observers to perform more efficient noise filtering at the price of an increase in controller design complexity. Remark 7. This paper does not address the measurement noise issue explicitly. The experience with the cascade MPC controller in [17] suggests that low and moderate measurement noises are not problematic, and simple low-pass filtering can be used. This is also favored in [33]. However, high measurement noises present a challenge to the implementation of linear static state-feedback controllers. The problem lies in the need for differentiation of the pendulum angle and the cart position w.r.t. time. In such a situation, it is useful to consider sophisticated differentiations or other controller structures such as the output-feedback controller with or without observers to perform more efficient noise filtering at the price of an increase in controller design complexity.

Conclusions
Most of the standard non-MPC linear static state-feedback controllers that are used to stabilize the nonlinear fourth-order CIP system around its desired equilibrium points are designed using pole-dependent techniques and without referring to a meaningful and transparent cost function. They are also prone to parameter tuning difficulties when dealing with the system transient performance and stability robustness simultaneously. To circumvent these difficulties, a linear static pole-independent MPC controller with only two tuning parameters is proposed. The controller ensures a priori a prescribed system gain margin and a two-time scale structure, allowing the response of the pendulum to be fast in comparison to the cart response. Using the linearized CIP model, some guidelines are developed in the form of two-dimensional contour plots for rapid MPC parameter tuning. Compared to the optimal CMD, CAS, and CPP methods, the proposed SIM method exhibits the best performance in the absence of disturbance input at the price of a slight degradation in the cart position overshoot and settling time. Future works will be devoted to studying the region of attraction associated with the proposed controller, reducing the sensitivity of the proposed method to disturbance input by optimizing the two-time-scale structure, relaxing the hypothesis of Assumption 1, and applying a similar method for other (linearized) robotic systems, e.g., the rotary inverted pendulum robot and the CIP with multiple links. The extension of the proposed method to deal directly with nonlinear system dynamics with or without uncertainties is also an interesting research topic. Solving (A5) gives the SFC gains as follows: where the GTC, i.e., τ, is the only parameter to be adjusted. On the other hand, the transfer function (A4), which has two real zeroes, can be decomposed as follows: where C 1 (s) is a non-all-pole system with one pair of jω-axis zeroes, and C 2 (s) is an all-pole system. Under the Manabe form, the authors of [28] showed the possibility to obtain a monotonic step response for C 1 (s) if τ is kept above a certain lower bound. However, obtaining the step response monotonicity for (A7) is almost impossible due to the presence of the NMP zero in this transfer function [29,30]. Now, combining (12), (22), and (A6) yields: where τ 0 = a 1/2 1 τ. The maximum system gain margin is equal to 0.5 + 2 √ 5/5 ≈ 1.39 at the normalized GTC τ 0 = 500 1/4 ≈ 4.73.

Appendix B. The Cascade MPC Controller (CAS)
In a previous paper [12], we have developed a two-horizon time cascade linear MPC controller to stabilize the CIP system. The controller scheme is composed of two MPC controllers in cascade. These controllers are tuned so as to obtain a doubly critically damped behavior for the inner and outer loops. For the purpose of a fair comparison, the cascade MPC (CAS) controller presented hereafter and shown in Figure A1 is an adapted version of the one developed in [12]. More precisely, to develop this version, we assume, as for the proposed SIMO MPC controller of this paper, that Assumption 1 is satisfied. In this case, the associated CAS pendulum prediction model, which is identical to the SIMO one, and the CAS cart prediction model, which is distinct from the SIMO one, are given hereafter: where (A i , B i , E i ) with i = 1, 2 are evaluated using (8), and x 1d is an unknown low-speed pendulum angle reference.  The cascade controller design can be conducted using two successive steps. In the first step, the design of the inner MPC controller ( )  The cascade controller design can be conducted using two successive steps. In the first step, the design of the inner MPC controller K 1 = (K 11 , K 12 , K 13 ) is considered using the following quadratic cost function: . (A11) From the angle prediction model (A9), the control input is obtained by minimization of (A11) as follows: v 0 (t) = −K 11 x 1 (t) − K 12 x 2 (t) + K 13 x 1d Imposing the critically damped constraint to the second-order inner loop gives [12]: where d 0 = 2a −1/2 1 ω n0 ≈ 1.2h −1 0 with the inner loop undamped natural frequency ω n0 ≈ 0.6h −1 .
In the second step, the design of the outer MPC controller K 2 = (K 21 , K 22 , K 23 ) is considered using a second quadratic cost function: J 2 (t, h 2 ) = e 2 3 (t + h 2 ) + ρ 2 e 2 4 (t + h 2 ) e 3 (t + h 2 ) = x 3d −x 3 (t + h 2 ) e 4 (t + h 2 ) = 0 −x 4 (t + h 2 ) . (A14) Using the cart prediction model (A10), the remaining control input is obtained by minimization of (A13) as follows: x 1d0 (t) = −K 21 x 3 (t) − K 22 x 4 (t) + K 23 x 3d Assuming that K 21 = K 22 = K 23 = K and ω n0 is high enough to obtain, the pole dominant principle leads to approximate the fourth-order cart transfer function by a simple second-order model [12]. Tuning the gain K to obtain a critically damped behavior for the outer loop leads to determining it with the following expression: From Figure A1, we get the following ECS: with where the cascade MPC controller gains are given in Equations (A13) and (A16). Finally, let us notice here that by construction, the obtained cascade MPC gains are positive. The trivial prediction models, defined in (A9) and (A10), appear to be sufficient to conduct the design, and there is no need to define a generalized prediction model as in the case of the SIMO MPC method. Additionally, the evaluation of the optimal system gain margin for this controller leads to solving a fourth-order polynomial equation that is preferable to solve numerically.

Appendix C. The Coincident Pole Placement Controller (CPP)
Taking into account the usefulness of driving the cart system with little or no oscillation to x 3d , a standard approach to design SFC controllers consist of determining the SFC controller gains (N 1 , N 2 , N 3 , N 4 ) so as to obtain a coincident real negative pole p structure for the targeted closed-loop system F 3 (s) [17]. In this situation, the CIP characteristic polynomial follows the configuration P(s) = (s − p) 4 that leads directly to evaluate the SFC controller gains by: The above controller can be viewed as optimal if its pole maximizes the gain margin (22). To this end, let us define the following normalized pole parameter: p 0 = a −1/2 1 p. (A20) Then, from (31), (A19), and (A20), we obtain the set of algebraic equations: Now, substituting (31) in (22) and combining the obtained result with (A21) yields: GM min = N 20 D −1 0 (r 12 − r 34 ) 1 + N 40 D −1 0 (r 12 − r 34 ) The maximum gain margin is equal to 0.5 1 + √ 5 ≈ 1.62 at the pole location p 0 = −5 −1/4 ≈ −0.6687. As it is stated in [17], the parameters (A19) that result from the above pole location choice allow one to tune the controller gains for any CIP system in the simplest fashion to ensure significantly improved performance and robustness. In this paper, we take this controller as a reference for the comparison task.